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Observation of epidermal cells or cell fragments on flat adhesive substrates 
has revealed two distinct morphological and functional states: A non-migrating 
symmetric "unpolarized state" and a migrating asymmetric "polarized state". 
They are characterized by different spatial distributions and dynamics of impor- 
tant molecular components as F-actin and myosin-II within the cytoplasm, and 
integrin receptors at the plasma membrane contacting the substratum, thereby 
inducing so-called focal adhesion complexes. So far, mathematical models have 
reduced this phenomenon either to gradients in regulatory control and signaling 
molecules or to different mechanics of the polymerizing and contracting actin 
filament system in different regions of the cell edge. 

Here we offer an alternative self-organizational model in order to reproduce 
and explain the emergence of both functional states for a certain range of dy- 
namical and kinetic model parameters. We apply an extended version of a 
two-phase, highly viscous cytoplasmic flow model with variable force balance 
equations at the moving edge, coupled to a 4-state reaction-diffusion-transport 
system for the bound and unbound integrin receptors beneath the spreading 
cell or cell fragment. In particular, we use simulations of a simplified 1-D model 
for a cell fragment of fixed extension to demonstrate characteristic changes in 
the concentration profiles for actin, myosin and doubly bound integrin, as they 
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occur during transition from the symmetric stationary state to the polarized 
migrating state. In the latter case the substratum experiences a low magnitude 
"pulling force" within the larger front region, and an opposing high magnitude 
"disruptive force" at the shorter rear region. Moreover, simulations of the cor- 
responding 2-D model with free boundary show characteristic undulating pro- 
trusions and retractions of the cell (fragment) edge, with local accumulation of 
doubly bound adhesion receptors behind it, combined with a modulated retro- 
grade F-actin flow. Finally, for a stationary model cell (fragment) of symmetric 
round shape, larger fluctuations in the circumferential protrusion activity and 
adhesion kinetics can break the radial symmetry and induce a gradual polar- 
ization of shape and concentration profiles, leading to continuous migration in 
direction of the leading front. 

The aim of the chapter is to show how relatively simple laws for the small- 
scale mechanics and kinetics of participating molecules, responsible for the en- 
ergy consuming steps as filament polymerization, pushing and sliding, binding 
and pulling on adhesion sites, can be combined into a non-linearly coupled sys- 
tem of hyperbolic, parabolic and elliptic differential equations that reproduce 
the emergent behavior of polarization and migration on the large-scale cell level. 
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1 Biology of cell polarization and migration 

Cell polarization and migration plays a central role in the development and maintenance 
of tissues in multicellular organisms. During ontogenesis, new tissues are formed by the 
coordinated division and locomotion of single cells. The polarization of a cell not only 
defines the direction of migration [3U] but also the cell division axis [M] and thus the 
three-dimensional structure of tissues, organs and finally the whole organism. 

1.1 Asymmetry of actin polymerization and substrate adhesion 

The coordinated development and release of focal adhesions (FAs) is a basic requirement for 
directed cell migration. Migrating cells feature pronounced adhesion dynamics and a struc- 
tural polarity with a clearly distinguishable frontal and rear area. Actin polymerization 
predominates at the cell front resulting in a protruding lamella in direction of migration. 
New focal adhesions composed of clustered protein complexes develop at the lower mem- 
brane of the lamella near the leading edge and couple the F-actin-network mechanically 
to extracellular matrix proteins. Simultaneously, the matured FAs residing at the opposed 
trailing edge are dissolved while myosin driven contraction of F-actin moves the cell body 
forward HQIES]. 

Arising from the process described above, the polarity of the cell can be referred to two 
structural asymmetries, which are key requirements for effective cell migration: asymmetry 
of actin polymerization and asymmetry of adhesion strength. The growth of actin filaments 
has to predominate at the cell front for pushing the leading edge in direction of migration 
|63l I58| . To move the cell body forward, it has to be released from the substrate during 
contraction by dissolving the FAs at the back (rear release) while the FAs at the cell front 
have to remain stable to provide a mechanical attachment for the contractile machinery 
pulling the cell body. In absence of rear release, traction forces could be dominated by 
adhesion forces and the cell would get stuck |54J. In this regard, the spatial distribution 
of adhesion strength and actin polymerization defines the direction of migration and could 
be specifically regulated due to directed cell movement. 
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Figure 1: Cell fragment experiments where the mechanical stress (release) induces a transi- 
tion from a circular non- migrating state (a) to a polarized migrating state (b,c). 
The cell fragment was labelled with fluorescent myosin II. Time scale is min:sec 
and length scale of bar is 2 fim (from |63j). 



1.2 Flow of actin filaments and myosin gradient 

The mechanisms underlying the structural polarity of migrating cells are still under discus- 
sion, particularly what concerns an intrinsic directionality of the cytoplasm. For example, 
the finding that asymmetric adhesive structures define polarization of a touching cell [30] 
suggests that directional flow of the actin cytoskeleton is involved in the process of cell 
polarization. Moreover, recent studies on fish keratocytes have revealed that polarization 
occurs spontaneously and is accompanied by a reorganization of the actin cytoskeleton, 
which finally leads to cell locomotion [64 . In these experiments, unpolarized solitary cells 
feature a circular shape with a radially symmetric actin distribution. Although the cells 
do not move, transient protrusions and retractions appear at the cell edge, and actin flows 
centripetally with a decreasing flow gradient from the cell edge to the center. In this ap- 
parently unstable state spontaneous symmetry breaking results in a faster inward flow and 
in an increased concentration of actin at the rear region of the cell. On the opposed front 
edge, the reduced inward actin flow causes protrusion and the development of a lamel- 
lipodium. In this polarized state, the cell starts to migrate while attaining a more or less 
constant shape. 

A similar behavior was observed with cell fragments extracted from the lamella of fish 
keratocytes [63], see Fig. [T] Since these fragments were lacking most cell organelles and 
the microtubule system, they were mainly limited to the actin-myosin machinery. They 
appeared either in an unpolarized, non-migrating state with homogeneous actin-myosin 
distribution or a polarized, migrating state with a rising myosin concentration from front 
to rear. Interestingly, polarization of fragments could be induced by mechanical stimula- 
tion, e.g. by shear flow or stress release leading to a transition from the stationary to the 
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migrating state. 

Cell polarization is discussed to result from the local activation of GTPases by the rho- 
family \62\ [54"1 143] . These proteins are known to regulate actin polymerization, myosin 
activity and FA assembly, whereas it remains unclear how the distribution of the GTPases 
is controlled. Recent experiments have revealed that by inhibition of myosin activating 
signalling pathways the cell's ability to polarize is reduced [53]. However, the observations 
of cells or cell fragments to polarize either spontaneously or due to mechanical pertur- 
bations, raise doubts that biochemical signalling is the primary reason for inducing and 
maintaining polarity. Experiments focussing on cell mechanics rather show that direction 
and strength of locomotion forces are inherently connected with the retrograde F-actin flow 
|53|. I22J . This suggests that the mechanical action of myosin, namely to induce cytoplasmic 
contraction, flow and force transduction at adhesive sites, plays a key role in explaining 
the ubiquitously observed phenomena of cell polarization and locomotion. 

2 Previous models of cytoplasm and adhesion dynamics 

The first detailed mathematical model coupling cytoskeleton and adhesion dynamics was 
developed by Lauffenburger and coworkers almost two decades ago [T7]. Here, the con- 
tractile actin-myosin network and its interaction with bound adhesion sites is represented 
as a mechanical system of connected viscoelastic units of generalized Kelvin- Voigt type, 
which constitute the (three) inner segments as well as a front segment (lamellipod) and a 
tail segment (uropod) for a rectangular model cell of fixed width and length. A coupled 
system of reaction-diffusion equations for free and substrate bound integrins or adhesion 
receptors on the dorsal (lower) and ventral (upper) part of the model cell is solved un- 
der pseudo steady-state assumptions. This dynamically provides the number of adhesion 
bonds in both end segments, lamellipod and uropod, whereby their rupture (dissociation) 
kinetics exponentially increases with force load onto a bound adhesion site. The resulting 
non-linear dynamics for the local (forward or backward) displacement of each viscoelastic 
unit produces a persistent forward translocation of the whole model cell. However, this 
is achieved only if a front-tail asymmetry is presupposed, either by exposing more free 
adhesion receptors, or by assuming higher adhesion bond affinity at the front compared 
to the tail. The authors present a series of results on how the simulated cell migration 
speed depends on various model parameters as, for example, cytoskeletal contractility or 
adhesion strength. A later variant of this model presents a more explicit study of adhesion 
bond disruption kinetics at the rear of the cell and already uses a four-state model for 
integrin binding to the cytoskeleton and to the substratum [51] . 

In recent years a series of more elaborate models have been developed, accounting for 
details of the meanwhile discovered molecular regulation mechanisms for the chemical and 
mechanical processes, particularly at the free boundary of a moving cell. One model type 
is based on spatially discrete algorithms (cellular Pott's model) using the definition of local 



5 



energies to determine the protrusion and retraction of boundary elements via the stochastic 
Metropolis rule, e.g. [331 lift] 127] . 

Another class of mechanical cell models takes into account branching and anisotropy of 
cytoskeletal actin filaments at various times and in various regions of the cell, see e.g. [65 [ I23 [ 
[28] . Besides explicit constructions of an elastically cross- linked network in the leading edge 
[33 [32] biophysical models have been developed to capture different dynamics of Arp2/3- 
induced branching and myosin-induced contraction of F-actin networks by Mogilner and 
coworkers [44 \ [25], [53] . The last article is particularly devoted to explain the mentioned 
polarization experiments with cell fragments (see Fig. [TJ by assuming a different actin- 
myosin network organization at the rear in comparison to the front region. 

Based on the 'reactive flow' model by Dembo and coworkers [T3l [12], [3] , the most elabo- 
rate model extension has been presented by Oliver et al. [49J: they use full 3-D two-phase 
flow equations with free ventral boundary and with two additional rapidly diffusing mes- 
senger concentrations that regulate actin network contractility and (de-)polymerization. 
Moreover, cell adhesion is modeled by a Navier-slip boundary condition at the substra- 
tum, in which only constant adhesive properties are taken into account. The analysis of 
this complex model, being performed in the thin film limit, is restricted to linear-stability 
arguments for ruffle generation and to local expansion analyses at the moving tip; there, 
phenomenological equations for boundary mass fluxes are considered without specifying the 
types of molecular mechanisms for tip protrusion. Finally, quantitative estimations for the 
pseudopod protrusion and cell translocation speed under various sub-limit assumptions are 
given, which turn out to be consistent with observed values, particularly for osteoblasts, 
though no numerical simulations are given that would re-insure the analytic results. A 
similar thin-film approximation of the 2-D equations under incompressibility assumptions 
for a 'viscous polar gel' has been used in [33] to derive explicit expressions for the advancing 
speed of a cell lamella, but again by predefining its polarity. 

Except the last one, all so-far mentioned models do not explicitly quantify the varying 
force field, which is applied by the cytoskeleton onto the substrate covered by a migrating 
cell and which has been approximately reconstructed by different inverse methods from 
experimental assays of cells, e.g. moving on flexible substrata [TH] [3H] [33 E|- A first 
model implementing force transduction to the substrate has been proposed by Gracheva 
and Othmer [24j by specifying a spatially one-dimensional system of viscoelastic equations 
for cytoskeleton dynamics, whose polymerization, contractility and adhesive binding is reg- 
ulated by signalling molecules. However, they make a pseudo-steady-state assumption for 
the binding kinetics of myosin polymers to actin filaments, and of transmembrane integrin 
proteins to the substratum. Moreover, they impose artificially defined gradients from tail 
to front of certain regulatory proteins in order to stimulate polarized cell translocation. 

Recently, adhesion kinetics has also been implemented into an extended cytoplasm flow 
model [31 E] describing the F-actin dynamics in an annular domain and its coupling to lamel- 
lipodial protrusions and retractions [59] . Force dependent maturation of FA complexes and 
active polymerization of F-actin enable the simulation model to reproduce characteristics 
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of fibroblast shape deformation and translocation. 

In an earlier publication we have presented another extension of the basic two-phase flow 
model for the 2-dimensional viscous cytoplasm dynamics [2j by coupling the constitutive 
hyperbolic-elliptic equations to a system of four reaction-diffusion-transport equations for 
the integrins beneath the cell or cell fragment [35]. Here we propose a generalized contin- 
uum model, in which we couple cytoplasm and adhesion dynamics with mechanical tension 
and transport of the plasma membrane, in order to reproduce spontaneous and induced 
cell polarization leading to migration, with assembly of adhesion sites at the cell front and 
adhesion release at the cell's rear end. Moreover, the model exhibits typical features of 
migrating cells as protrusion/retraction cycles, rearward actin flow, pulling forces at the 
front and a concentration of disruptive forces at the rear. In the model, the cytoplasm is 
described as a viscous and contractile fluid of polymers representing the actin cytoskele- 
ton interpenetrated by an aqueous phase. This actin filament network is preferentially 
assembled at the cell edge, and can be contracted by cross-linking with diffusing myosin 
oligomers. The moving actin filaments then couple to transmembrane adhesion proteins 
that are freely diffusing in the membrane or bound to the substrate on the extracellular 
domain. Thus, the cytoskeleton mechanically connects to the substrate through dynamic 
binding processes which results in force transduction and finally cell locomotion. 

Throughout our model presentation we rely on continuum descriptions, in which macro- 
scopic mass and momentum laws are combined with mesoscopic submodels for fast molec- 
ular kinetics due to adequate pseudo-steady-state assumptions. 

3 Two-phase flow model for cytoplasm coupled to reaction, 
diffusion and transport for myosin-ll and integrin proteins 

For simplicity, we restrict our model derivation and analysis to a flat two-dimensional ge- 
ometry, so that cells or cell fragments are assumed to be homogeneously spread on the 
substrate without considerable change of cell height. Thereby, three-dimensional effects 
around the cell nucleus (e.g. due to cell rolling) or along the ventral (upper) plasma mem- 
brane are neglected. In order to reproduce the main biophysical mechanisms and biochem- 
ical processes that enable a cell to polarize and translocate on an adhesive substratum, we 
nevertheless distinguish between the cytoplasm and the exterior plasma membrane: On the 
one hand, F-actin assembly, myosin kinetics and viscous network flow take place within 
the cytoplasmic interior of the cell. On the other hand, integrin binding to substrate 
or cytoskeleton and its transport or diffusion are confined to the dorsal (lower) plasma 
membrane, where forces are transduced to or generated at the cell periphery, leading to 
tip protrusion or retraction. Thus, the moving cell (fragment) is simply represented by 
a time-dependent connected domain Q(t) C M 2 , over which the cytoplasmic volume ex- 
tends with fixed constant height, and any ruffles or blebs on top of the cell are neglected. 
We rather assume that the cell dynamics is completely determined by a flat cytoskeleton 
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Figure 2: Schematic longitudinal section through a cell fragment (or a cell lamella) of con- 
stant height with involved proteins. For further explanation and model derivation 
see the following sections. 



sheet of F-actin network, see Fig. [2j This cross-linked filament phase with volume fraction 
6(t, x),x £ Q(t) is of constant thickness and connected to the upper and lower plasma 
membrane in such a way that also the cytosol, i.e. the solvent phase with volume fraction 
(1 — 9(t, x)), is confined to the same volume space. 

By suitable model simplifications, we aim to capture the self-organizational power of the 
cytoplasm as a two-phase fluid coupled with the reaction, transport and diffusion of a series 
of chemical ingredients. We explicitly model the kinetics and dynamics of only those F-actin 
associated proteins that induce the main bio-mechanical processes of force generation and 
transduction to the substratum, namely myosin oligomers and transmembrane integrin 
proteins. All other regulatory proteins such as e.g. Rac, Rho or the branching protein 
Arp2/3 or smaller substrate molecules, e.g. monomeric G-actin, are assumed to rapidly 
diffuse within the cytosol under fast regeneration, so that they attain constant reservoir 
concentrations serving as parameters in the model. 

3.1 Mass balance and flow equations 

In the chosen two-dimensional continuum description, cell deformation and transloca- 
tion with respect to the fixed substratum is represented by movement of the cell edge 
r(t) = dQ(t) which, however, is induced by three distinguished mass flows with poten- 
tially different transport velocity fields on 0, namely the mean velocities v of the F-actin 
cytoskeleton (6) and w of the aqueous cytosol (1 — 0) plus the transport velocity u of the 
lower dorsal plasma membrane. Thereby both v and w are averaged over the constant 
cell height. For the membrane area flux u we make the simplifying assumption that the 
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lipid bilayer between the flat substratum and an adhering cell is always stretched out to 
its maximal extension with constant area concentration of the lipid-protein mixture, thus 
constituting an incompressible two-dimensional Newtonian fluid. Then, freely diffusing 
transmembrane integrin proteins are additionally transported by the membrane velocity u, 
myosin oligomers by the cytosol velocity w, and F-actin bound myosin-II motor molecules 
by the cytoskeleton velocity v. The detailed mass balance equations are discussed in the 
following paragraphs. 

3.1.1 Mass conservation for the two phases cytoskeleton and cytosol. 

Cytoskeleton and cytosol in the flat two-dimensional geometry can experience counteracting 
flows due to local contraction, assembly or relaxation of the F-actin network, while the bulk 
cytoplasmic fluid with constant volume fraction 1 can be assumed to be incompressible, 
at least in the range of occurring pressure differences < kPa. Compare the analogous 
situation of a water-filled sponge that is internally condensing without changing its shape. 
Then, because of the fixed height assumption, the total two-dimensional volume flux W = 
#v + (1 — 0)w is divergence free at any time t, yielding the first local mass conservation 
law 

V- (0v + (l-6)w) =0. (1) 

Therefore, also the total cell volume (i.e. the 2-dimensional cell area) is conserved over 
time, so that the 'bulk' fluid moves together with the free boundary meaning that 

the total volume flux W has to fulfill the natural free boundary condition 

u v . W = f (2) 

with vy denoting the exterior normal of T and T quantifying the normal speed of the cell 
edge. 

In addition to possible convection, the F-actin network can locally be assembled by 
filament polymerization from the pool of monomeric G-actin within the cytosol and it 
can be disassembled by the reverse process of severing or depolymerization. This mass 
exchange between the two phases is given by a second conservation law, which for the 
filament volume fraction 6 can be written as 

dtO + V • (0v) = R{9) (3) 

with a net assembly rate R{6) to be modeled. According to our simplifying assumption, the 
dependence on G-actin, Arp2/3 and possible regulatory proteins enters only via constant 
parameters: 

R{6) = {k on a g - fc off ) B(9) - r6 (4) 

with concentration of globular actin a g , polymerization rate k on and depolymerization rate 
fc ff & t the filaments' plus (barbed) ends, as well as another lumped disassembly rate r, 
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see e.g. [12]. Due to relatively fast nucleation and capping of actin filaments, the relative 
number B of barbed ends is assumed to stay in pseudo-equilibrium with the local F-actin 
concentration a = 6 ■ a maX ; where we suppose a maximal condensation of actin filaments in 
the order of a max = 800 /jM. Following [32] we write 

B(6) = ^(e + v Q ° ) (5) 

with capping rate u, a basic nucleation rate e and an induced branching rate v§ = f a Arp0, 
proportional to the concentration ArpO \pM] of activated Arp2/3, together with a half- 
saturation concentration K a for its primary actin binding site. 

3.1.2 Reaction-transport-diffusion equations for myosin oligomers. 

Myosin-II oligomers are the most important actin binding proteins that are responsible for 
the generation of contractile forces within cross-linked F-actin networks. Thus, their spatial 
distribution within a polarizing or moving cell plays a key role for the cytoplasm dynamics. 
In a most simple way we only distinguish between freely diffusing myosin oligomers (nif) 
and those that are bound to cytoskeletal actin filaments (mt) and, therefore, are convected 
with velocity v. Since free myosin tetramers first have to attach to the actin network at one 
binding site, we consider this bimolecular reaction (m/ with a) as the rate limiting step. 
Correspondingly we assume, that the relatively faster processes of double binding (m b with 
a) and power stroke formation are in a pseudo-steady state. Thereby the contractile stress 



ij) = with 9 = a/a max in the network is determined, see equation (25 ) in the following 

section. Finally, we suppose a constant diffusivity D m for free myosin oligomers, embedded 
into the cytosol flow w. Under these assumptions the local mass balance equations for the 
corresponding myosin concentrations are: 



d t mf = V • [p m Vmf - mjwj - a m ■ a ■ m/ + 5 m (a) ■ m b (6) 

dtnib = -V • (m b v) + a m ■ a ■ m f - 5 m (a) ■ m b . (7) 

For the association rate a m we assume a constant parameter, whereas the dissociation rate 
is supposed to increase quadratically with increasing network concentration a because of 
steric inhibition or competition for binding sites: 

M<0 = &(i + ^r) =^( 1 + ^;)' ( § ) 

with an optimal actin network concentration a op t = Omax^opt- Moreover, in the case of fast 
diffusion and no transport, the constant equilibrium concentrations m*^ and satisfy the 
equality 

* OL m a ^ . . 

o m [a) 1 

so that the resulting contractile stress ij)(0^) becomes optimal for 9 = $ ptj 

Fig. [5] below. 
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3.1.3 Mass conserving flow of the dorsal plasma membrane. 

As mentioned before, we will assume that the dorsal plasma membrane beneath an adhering 
cell (fragment) is stretched and that small fluctuations can be neglected. On the contrary 
top side, the ventral plasma-membrane usually shows extensive wrinkles or folds, which 
most probably are induced by the contractile action of the cytoskeleton itself, indicating 
that the plasma membrane tip at the cell edge stays under a positive tension T r . Before 
discussing the corresponding dynamics of the dorsal membrane, we state the two extreme 
possibilities, namely whether the membrane moves together with the cell or not: 

1. Membrane sticking to substrate: During cell translocation over the substratum the 
whole dorsal membrane stays fixed and no slip may occur due to relatively strong 
interaction forces with the substratum, e.g. via the glykocalix, so that u = 0. 



2. Membrane sticking to cell edges: The dorsal membrane is pulled along the substratum 
due to cytoplasm protrusions at the leading front, but with no slip around the lamellar 
tips due to strong membrane curvature. In this way the membrane area flux satisfies 
the mass conservation law together with a boundary condition analogous to (|2j), 
namely 

V-u = onO(t) v T ■ u = f on T(t). (10) 

In the latter case, the membrane can slip over the substratum, thereby experiencing a 



finite frictional drag force, see (40) below, whereas in the first case this force is infinitely 
large to suppress slipping - but then the membrane has to slip around the tips at moving 
cell edges. Clearly, a certain mixture of both possibilities is physically realizable but not 
considered here. 

Moreover, frictional drag onto the dorsal membrane can also occur at its cytoplasmic 
side. Whereas the relative flow of cytosol will have negligible effects, the horizontal F-actin 
flow develops a vertical flow profile. This profile depends on the amount of 'Navier-slip' 
in a cortical layer on top of the dorsal membrane and emerges due to frictional forces 
between actin filaments and membrane proteins. In analogy to the more explicit thin film 
approximation |49j we only consider the averaged profile velocity v. In this way we can 
introduce a simplifying cortical slip parameter, < k < 1, so that the effective relative 
velocity between cytoskeleton and dorsal membrane is reduced to k(v — u). Thus, the 
resulting effective velocity of cortical F-actin in a thin layer above dorsal membrane and 
substratum is given by 

v c = kv + (1 — ft)u. (11) 

Here the factor k, which clearly depends on the viscous shear properties of the cytoskeleton, 
should be larger than 0.5 to reflect observations of a generally slippy behavior (see e.g. [29^ 
I22j). Moreover, we suppose that the vertical profile, thus v c , is not changed if some of the 
cortical actin filaments are (transiently) bound to integrins in the dorsal membrane, see 
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Figure 3: Scheme of four states of the transmembrane adhesion protein integrin. (/): un- 
bound and freely diffusing within the dorsal membrane, (s): bound to the sub- 
strate, (a): bound to the actin network and thus moving with the actin-flow, (sa): 
bound to actin and substrate (force-transducing state). Integrins can switch be- 
tween these states due to reversible binding kinetics with binding/unbinding rates 
a, f3, 7 and 5. 



Fig. [3] and the following paragraph: Either those bound integrins or integrin complexes are 
passively pulled through the lipid-protein bilayer with relative velocity v c — u or, in case 
of substrate-fixed adhesion bonds, with relative velocity — u. The whole F-actin network 
above such a bond is slowed down by a certain local frictional force per adhesion site, 

F c = $ o 0v c , (12) 



entering into the corresponding force balance law, see (23) below. In any case, the fric- 
tional drag onto the dorsal membrane, induced by the relative motion of singular integrin 
complexes, will be neglected in our model. 



3.1.4 Reaction-transport-diffusion equations for membrane integrins. 

The mechanical connection between the actin cytoskeleton and extracellular matrix pro- 
teins is provided by transmembrane integrins appearing in four different states [51j [35] , see 
Fig-i 

In the freely diffusing state /, integrins are neither coupled to the substrate nor to the 
actin cytoskeleton and move according to a simple diffusion-transport law within the dorsal 
membrane. These integrins can change their state by binding to the actin cytoskeleton (a) 
or to the substrate (s). In state a, the integrins move with the cortical F-actin velocity 
v c , whereas integrins in state s remain stationary with respect to the substrate. Actin- or 
substrate-bound integrins can switch back into the freely diffusing state / by unbinding, 
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or into the state sa by coupling to the cytoskeleton and the extracellular matrix. Cell ad- 
hesion occurs only in this double bound state s representing a focal adhesion (FA), where 



the frictional force F c (12) of the moving actin-network is transduced to the substrate. 
The concentration of integrins c# in the different states is described by the following cou- 
pled system of differential equations consisting of terms for spatial movement and binding 
kinetics. 

d t cf = V • (DfVcf - c/u) + a~c s + 5~c a - (S^a + a + )c/ (13) 
d t c a = -V • (c a v c ) + S^acj - (5~ + 7 + )c a + 7^ exp(p 7 |F c |)c sa (14) 
d t c s = a + Cf - (a~ + f3^a)c s + /3q exp(p /3 \F c \)c sa (15) 

d t c sa = 1 + c a + PqClc s - exp^lFd) +7^ exp(p 7 |F c |)^c sa (16) 

The bonds of the force-transducing integ rins (c 5a ) can break when they experience a me- 
chanical stress by the cytoskeleton, in our model given by the modulus of the force vector 



F c defined in (12). Then, according to the theory of Bell [57] the dissociation rates 7 



and j3~ depend exponentially on the mechanical load |F C |, see equations (14-16) above, 
with 7q~ and (3q describing the basic dissociation rates without load and the exponential 
coefficients p# = C#/ksT measuring potentially different rupture rates from the substrate 
or cytoskeleton binding site, respectively. 

3.1.5 Mass flux conditions at the free boundary. 

In addition to the already mentioned bulk flux conditions ^ and (11) related to the 
normal speed of the cell edge T(t), we have to impose compatible boundary conditions 
onto the concentrations of those molecular species that are not fixed to the substratum. 
For the two parabolic diffusion equations ^ and (13) we impose, respectively, natural zero- 



flux boundary conditions onto freely diffusing myosin (m/) and fixed Dirichlet conditions 
(cf = Cj) onto freely diffusing integrin. Thereby we suppose a constant reservoir of fresh 
adhesion receptors expressed in the upper ventral membrane and diffusing (or eventually 
being transported) from there around the lamellar tip into the lower dorsal part of the 
membrane. 

The F-actin flux 9v cannot leave the cell, meaning that on the free boundary T the 
relative inward normal F-actin velocity always has to satisfy the inequality 

F = r- i/r-v>0. (17) 

If the strong inequality V > holds at a certain boundary point, then two different 
modeling situations may arise for the cell edge: 

1. No-stick condition at the lamellar tip: Local disruption of the actin network from 
the edge is allowed, e.g. under suitable load conditions [35], so there is no new F- 
actin production directly at plasma membrane edge. Therefore we have to impose 
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zero-Dirichlet condition for the F-actin concentration 



9 = if V > holds at a non-sticky point x S T. 



(18) 



Network sticking to the lamellar tip: As indicated in Fig. [4| active polymerization of 
actin filaments directly at the plasma membrane is allowed either (a) at fluctuating 
free filament ends touching the tip membrane [43J, or (b) at filaments that are bound 
to clamp- motor proteins anchored in the tip membrane [IB]. Both cases could occur 
simultaneously in a local region, but whether active polymerization with inward mass 



flux 9V > can take place depends on local force balance conditions, see Section 3.2 



and (36) below. 



Actin polymerization 

clamp-motor (a' : ) brownian ratchet (a } ) bound membrane 

protein (a h ) 




myosin 



actin 
monomer 



Figure 4: Schematic top view onto the lamellar tip of a cell (fragment) showing three pos- 
sibilities for actin filament ends to interact with the free plasma membrane: an- 
choring at a membrane protein, as integrin; freely fluctuating with polymerization 
by a "Brownian ratchet" mechanism; binding to a "clamp-motor" protein, like 
the WASP-complex, with induced polymerization. The resulting local pseudo- 
equilibrium F-actin concentrations are a b ,a^ , and a c , respectively. See equations 
(31) and (33) for a quantitative model. 



Finally, also the two hyperbolic equations ([7]) and ( 14 ) require zero-influx conditions, so 
that for the transported concentrations nib and c a we have to impose zero-Dirichlet condi- 
tions only in cases of V > or V c > 0, respectively. In the latter case the inward normal 



velocity of v c is defined in analogy to (17), noting that under the modeling hypothesis (10) 
we anyway have V c = kV. 
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3.2 Force balance equations 

3.2.1 Two-phase flow equations for cytoskeleton and cytosol. 

Since cell movement usually occurs on a time scale of minutes, and cell sizes usually are 
in the order of tens of microns, already the cytosol, with its consistency similar to an 
aqueous viscous fluid, has low Reynolds number. The more this is true for the cytoskeleton 
with consistency similar to a dense viscoelastic gel. The viscoelastic properties of the 
cytoplasm have been experimentally studied [X8|. [T9] and also mathematically modeled 
I49j . whereby the simple model of a Maxwell fluid seems to serve as a quite good 
description. This means, that elastic properties are mainly effective on a shorter time scale 
of seconds and dominated by viscous effects on a medium scale of minutes. Moreover, in 
contrast to a passive elastic material, the cytoplasm mostly stays under active contractile 
stress as exerted by myosin-II oligomers (in the order of kPa) , which is able to break weaker 
cross-links of the F-actin network e.g. by a-actinin or filamin, so that the contraction forces 
are equilibrated only by viscous and frictional forces. 

Under these assumptions, the highly viscous two-phase creeping flow model for cyto- 
plasm, originally proposed 25 years ago by Dembo et al. jE] and so far extensively ap- 
plied, see e.g. [U [26] , can be condensed into a pseudo-stationary linear elliptic system for 
the mean F-actin velocity v and the effective hydrostatic pressure p on O(t): 



V • /ztfVv + V ( S(9, m b ) -pj = $(0, c sa ) v c , (19) 

Vp = Lp6(v - w) . (20) 



Here the generalized (elliptic) Stokes equation (19) involves the effective stress S(9,nib 



as defined in equation (24) and the symmetrized displacement rate Vv (see e.g. [35J 



Moreover, v c is defined in (11), and /U,3> denote the coefficients of bulk viscosity and 
substratum friction, see (23). Because of negligible cytosol viscosity the simple Darcy law 
(20) suffices as a model description, where the coefficient <p measures the internal two- 
phase flow friction. From the last equation one explicitly solves for the cytosol velocity 
w = v — (l/ip0)S7p, so that the total volume flux is 

W = v + (l-0)(w- v) = v- ^— — Vp (21) 

Insertion into the mass-balance equation ([!]) then yields the generalized (elliptic) Laplace 
equation 

V • ^—^-Vp = V • v. (22) 

Notice that here the ellipticity degenerates for marginal volume fractions < 9 < 1, which 
can be relevant in cases where — > at boundary points of cytoskeleton disruption, see 
(18) above. 
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We remark that together with the hyperbolic mass balance equation ([3| the linear elliptic 
system ( 19|22 ) constitutes generalized pseudo-stationary Navier-Stokes equations for the 



F-actin network as a compressible, highly viscous and reactive fluid. It is mutually coupled 
to the mass concentrations in equations (I7|and (16) via a myosin-mediated contractile 



stress term appearing in S(9,rrib), see eq. ( p4h below, and an adhesion-mediated friction 
coefficient 

$(9,c sa ) = <$> c sa 9. (23) 



In this way the right hand side of (19) reads F„ = <i>v c = c sa F c , with the frictional force 



F c per doubly bound integrin defined in (12) above 



Tracing these model equations back to their derivation |X2(, [2l [35] provides us not only 
with precise biophysical conditions on v and p at the moving boundary T (see next para- 
graph), but also with genuine nonlinear parameter functions being based on thermody- 



namical reasoning at the molecular scale: The function S in (19) represents the effective 
stress in the network phase, which is induced by molecular interactions between the F- 
actin filaments themselves as well as by thermic interactions with solvent molecules in the 
cytosol. It is generally expressed as the weighted negative sum S = —9P a — (1 — 9)P S of 
corresponding pressures P a and P s which are applied to any network volume element and 
any cytosol element, respectively. Passive elastic stresses may be neglected since these are 
already relaxed on the creeping flow timescale. When finding a free binding site on fila- 
ments, previously single bound myosin-II tetramers exert power stroke motor forces with 
their free myosin heads, cf. Figs. [2] and |4} Thus, an attractive stress — P a = V'o'^b is applied, 
where the simple coefficient ipo comprises binding affinity, power stroke probability and the 
mean applied force per stroke. On the other hand, standard Gibbs free energy arguments 
suggest a molecular solvent pressure P s = — <roln(l — 9)/{l — 9), see p.]. Then we arrive at 
expressions for the 

effective stress S = S(9,mb) = tp(9,mb) — a(9) (24) 

contractile stress ip(9, mb) = ipo9rrib (25) 

swelling pressure a (9) = ao\ ln(l — 9)\, (26) 

see Fig. [5] for corresponding plots with nib = m % ™ equation ([9]). We emphasize that the 



effective stress function in equation (24) sums up the contributions from the cytoskeleton 
and the cytosol, so that both phases are not separated anymore. However, when deriving 
the corresponding natural boundary conditions for the elliptic problem, then the two phases 
have to be considered separately again. 

3.2.2 Stress and pressure balance conditions at the free boundary. 

Recall that according to ^ the cell edge T(t) always moves together with the normal 



component of the total cytoplasm flux W as quantified in (21). Then, from (17) also the 
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Figure 5: Plots of F-actin model functions. 

Dashed curves: net assembly 
rate R(9), equation Q, and ef- 
fective stress S(6) = tp(6) — 
cr(0), equation (24). Contin- 
uous curves: contractile stress 
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relative inward normal F-actin flux is determined as proportional to the inward normal 
pressure gradient on F(t) 

0V = e(t -z/r-v) = - 1 ^ d "r- Vp. (27) 

This explicit linear relation between free boundary speed, normal F-actin velocity and 
pressure gradient serves as an extra free boundary condition in addition to the set of 



boundary conditions that are necessary for uniquely solving the linear elliptic system (19) 



and (22) on the given domain Q(t). These depend on the particular model choice of 



boundary pressure functions at the cell edge membrane, but also on the outcome of the 



F-actin flow requirement V > in (17). 



Using the general derivations in [2J, eqs. (58-62), one obtains separate pressure balance 
conditions for each of the two phases at all free boundary points of T(t) satisfying V > 0, 

cytoskeleton —up ■ T Q ■ z/r + 6p =0P^ (28) 

and cytosol (1 - 6)p + a{6) =(1 - 6)P* (29) 

with the intrinsic stress tensor T a = //#Vv + ip(6, m&)I. These equations mean that at 
those parts of the tip plasma membrane, which are exposed to the cytoskeleton network 
(8) or to the cytosol (1 — 0), respectively, the sum of internal pressures is in balance with a 
certain boundary cytoskeleton pressure P^ or boundary cytosol pressure Pj at each volume 
element. Modeling expressions for these pressures will be given in the following paragraph. 



Summing up both pressure balance equations (28) and (29) yields the general Neumann- 
type condition for v on T(t) 

ur-Ta-ur- (a(8) + p) + 6P^ + (1 - = 0, (30) 
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which is now valid without restriction, see [2], eq. (58), i.e. also in boundary points where 



the condition V = holds. However, then the Dirichlet-type boundary condition (29) for p 



has to be replaced by a corresponding inequality, and insertion of V = into (27) provides 
a Neumann condition for p instead. In particular we have &r ■ Vp = 0, which also holds in 
case of disruption with V > and 9 = 0. 

Thus, in any constellation of the mentioned conditions on T(t), for fixed time t and given 
profile of all protein concentrations, we obtain a well-posed elliptic boundary value problem 
for v and p to be solved on (eventually by numeric iteration and/or conjugate gradient 
method, cf. |35j). Finally, the normal speed T of the free boundary can be calculated 
by using equation ([2Tb. This is also possible in the case of network disruption at some 
boundary point xo € T(t) with 0(t,xo) = 0, since then V(£, xo) > is obtained as the 
L'Hospital-limit of — 1 u r ■ Vp(£, x) as x — > xo from the interior. 

3.2.3 Boundary pressure functions at the cell edge. 

We only consider creeping cell migration in a medium that is not under external pressure or 
stress. Therefore, in a purely physical model, the boundary pressures and P s r along the 
cell edge F(t) should be set to zero. However, the edge is physically defined as the lamellar 
tip, where the ventral and dorsal part of the surrounding plasma membrane meet. Thus, 
if Tr denotes the scalar tension value of the dorsal plasma membrane at the cell edge (see 
the following paragraph) then this tension appears as a boundary tip pressure acting on 
both phases, the cytoskeleton and the cytosol. Moreover, in addition to the counteracting 



cytoskeleton and cytosol pressure at the left hand sides of eqs. (28) and (29), there could be 
extra pressures due to active polymerization forces, which usually are generated at barbed 
actin filament ends that can become exposed to the tip membrane in two variants, see 
Fig. [4] above: 

• Brownian ratchet model: Assume that at a point of the lamellar tip T(t) a frac- 
tion a B of filaments is bound to membrane proteins in a fast pseudo-steady state 
equilibrium with the actual F- actin concentration a = #a max . Then from the re- 
maining network with concentration a F (a) = a — a B (a), Arp2/3 induced branching 
can occur. The barbed filament ends are more or less normally exposed to the tip 
membrane with concentration a>f = a? (a) = «q • ArpO • K a +a F with suitably chosen 
coefficients, see also ([5]). Then, due to insertion of G-actin in-between fluctuating 
filament and membrane there appears a free polymerization pressure 



(31) 



with a ratchet coefficient tTq > depending on the free energy of one monomer 
addition, see relevant force estimations e.g. in |41j . 

Clamp-motor model: Alternatively or additionally, a certain fraction of tip bound 
actin filaments with concentration a c = a c (a) < a B {a) can be bound to WASP-like 
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membrane proteins, which serve as filament-end-tracking motors. By means of an en- 
ergy consuming polymer elongation process at the clamped end of the filament, these 
proteins push the bound filament outward, thus leading to a clamp polymerization 
pressure 

P c = tt c q - a c {a). (32) 

Here again, the clamp-motor coefficient 7Tq > depends on the energy of one monomer 
translocation, see [To] . 

The sum of these active polymerization pressures would induce an averaged relative normal 
inward mass flux OV > of the whole cytoskeleton, based on the particular inflows of a' - 
and a c -filaments satisfying 9V = a^Vf = a c V c . However, this flow experiences a viscous 
resistance from the remaining fixed filaments (concentration a B - a c ), with a viscosity 
coefficient that could increase in the presence of bound myosin molecules. Thus we get the 
total tip polymerization pressure 

-Ppoly =P f +P C ~ fJ-r (m b ) (a B - a c ) a max F. (33) 
Clearly, here the relative inward velocity V is the one defined by the normal component of 



v in relation to T, see (17). For a similar resistance model based on elastic cross-linking 
see [23]. 

Finally, since this intrinsic boundary pressure -P po i y acts onto the cytoskeleton volume 
fraction, while simultaneously inducing a corresponding counter-pressure onto the cytosol 
volume fraction, we can state the following distribution for the boundary pressures gener- 
ated at the membrane 

OPl = 9(1 - Kr )r r + Ppoiy , (34) 
(1 - 9)P? = {1-6 + 9K r )T r - P poly , (35) 

where we introduce a weight factor, < Kr < 1, measuring the relative effect of membrane 
tension rr onto the cytosol phase and possibly depending on the not-explicitly modeled tip 



geometry. Then, by substituting p from (29) into (30) we obtain the generalized Neumann- 
type boundary condition for v in all points of T(t) where the network is attached, a B > 0, 
and where V > holds: 

u T ■ T a ■ w - a{6) - ^Mm b ) (a B - a c ) V + (a{9) +p f + P c - 6k t t t ) = 0. (36) 

This means, that the mean inward F-actin polymerization speed V > on the moving cell 



edge r(t) is implicitly determined by solving the linear elliptic system (19) and (22) for 
v and p and satisfying all boundary conditions. Thereby the Neumann condition above 
contains all membrane protruding pressure terms in series, namely the swelling pressure a, 
the polymerization pressures pf and p c induced by a Brownian ratchet or an clamp-motor 
mechanism, as well as a counteracting stress due to dorsal membrane tension rr- 
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3.2.4 Global force balance at the adhesive substratum. 



There are two kinds of forces exerted by the migrating cell (fragment) onto the fixed flat 
substratum: the integrin-adhesion mediated active frictional force represented by the vector 



field F^ on the right hand side of Stokes' equation (U9L and an analogous passive frictional 



force F n due to motion of the dorsal plasma membrane along the substratum: 

F„ = c sa F c = $(e,c sa )v c , 
F u = $ n u, 



(37) 
(38) 



with <J> as in (23) and an assumed friction constant $ u > 0. Supposing that on the 



substratum no other forces are applied than these, then their integral sum has to vanish, 
so that the zero force balance holds: 







n(t) 



(F„ + F tt ). 



(39) 



Furthermore, if we assume that the dorsal membrane, viewed as a 2-dimensional incom- 



pressible fluid satisfying the zero-divergence condition (10), has relatively low viscosity. 



then the membrane tension t u induced by the frictional flow can be defined according to 
Darcy's law 

Vr u = $ u u (40) 

and thus determined as solution of the Laplace equation At u = with Neumann boundary 
condition up ■ Vt„ = 0, see (10). Then the zero force balance (39) together with (19) 



implies that the boundary tension values rp = r u |r(t), which are uniquely determined up to 
a constant, necessarily fulfill the integrability condition J r , t Jn6Vv+S(6, ^b)—p+ T T)^v = 0, 
which by insertion of the Neumann boundary condition ( |30| ) and by using the symmetry of 
Vv reduces to the equivalent necessary condition f r ^(0P£ + (l — 9)Pg —Tp)^ = 0. Indeed, 
this condition is fulfilled for the boundary pressure model functions that were chosen in 



the preceding paragraph, (34) and (35), because then even the integrand in the previous 
condition vanishes. 



4 Results of model simulations 

4.1 Spontaneous cell polarization in the 2-D model 

We have simulated the adhesive motion of a flat cell or cell fragment represented by a 
2-dimensional domain, £l(t), with moving cell edge or lamellar tip, T(t) = d£l(t), under 
certain simplifying assumptions: 

1. The lower dorsal membrane sticks to the substratum (u = 0) and there exists a small 
membrane tension rr > 0, constant over the whole cell edge. 
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2. There occurs no active polymerization pressure at the cell edge (pf = p c = 0), 
only the similar swelling pressure <j{6) and the hydrostatic pressure p can push the 
boundary. 

3. Disruption of the F-actin network from the lamellar tip T(t) can locally occur if the 
network tension exceeds a certain threshold that might depend on the fraction of 
membrane-bound actin filaments a B = A „ a . . 

K B +a 

4. The free myosin-II concentration is at a fixed constant level > and the amount of 
F-actin bound myosin-II oligomers is in a pseudo-steady equilibrium ml (a) according 
to ([9]), so that the contractile stress is only a function of = a/a max : 

m = Mmt = 5 o f l + 9 2 /6 2 • (41) 

More details and a list of chosen parameters can be found in |35j (Section 1 and Table 2). 
Since unpolarized cells and cell fragments, as observed under various conditions [63, 64J, 
attain a quite regular circular shape, we choose as initial condition a circle f2(0) of radius 
6 /mi. Moreover, in order to mimic the radial spreading of cells after exposure onto a 
flat substratum, we start with constant integrin densities and radially symmetric initial 
configuration for the volume fraction 9, with slightly larger values closer to the center. 
Due to this initial perturbation, the F-actin concentration rapidly condenses into a central 
region of high 6, surrounded by a lamella-kind region of low 8, see Fig. [6^a,b). Later, 
also the actin- and surface-bound integrin adhesion proteins c sa concentrate around this 
center, see Fig. [6^e). Both phenomena are supported by a strong retrograde F-actin flow, 
which collects actin filaments and actin bound integrins c a in radial direction from almost 
everywhere in the periphery. Moreover, the hydrodynamic pressure has its maximum 
within the center region (data not shown), so that its negative outward gradient represents 
the squeezed flow of cytosol out of the contracting F-actin network. 

At the free cell edge with relatively low F-actin concentration, there occurs a meta-stable 
equilibration between a positive swelling pressure pushing the lamellar tip outwards and 
a resisting viscous network tension pulling the tip inwards. For some time, see Fig. [6^b) 
after 5 minutes, local regions of protrusion or retraction can be observed, which point 
into varying directions along the cell periphery. Notice that these spatiotemporal fluctu- 
ations are not due to the tiny stochastic perturbations that are imposed to the F-actin 
polymerization rate, but represent the emergent chaotic dynamics of the cytoplasm as a 
reactive and contractile two-phase fluid [11] 15]. Furthermore, behind a cyclically protrud- 
ing and retracting free edge we can observe a layer with slightly increased concentration 
of substrate-and-actin bound integrin c sa , see Fig. |6je) . This is the onset of polarization: 
fresh free integrin proteins are appearing at the protruding part of the edge from the upper 
membrane, thus also increasing c sa in this region. As a consequence, a larger frictional 
force of the retrograde F-actin is generated in that direction, inducing a bias into the force 
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Figure 6: 2-D simulation of a spontaneously polarizing cell (fragment) starting with radially 
symmetric initial conditions: Spatial distributions shown as pseudocolor plots of 
F-actin volume fraction 9 in (a,b,c), and concentration of actin-and-surface bound 
integ rin proteins c sa in (d,e,f) at three different time instants: (a,d) 1 min, (b,e) 
5 min, and (c,f) 10 min after initialization. The width of the shown region is 
22 fim. 
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vector field transduced to the substratum. Thus, the whole cell fragment starts to move 
in this direction, which later becomes the leading edge of the migrating cell fragment, 
see Fig. J6jc,f). The emerging polarization of the cell is most clearly expressed in the c sa 
distribution of Fig. [6^f), where an increasingly dense band is formed behind the leading 
edge and the central region of focal adhesions is slowly shifted rearwards with increasing 
migration speed, while becoming deformed in shape similarly to the whole cell fragment. 

This simulation is an example for the autonomous formation of a meta-stable unpolar- 
ized, almost circular state and its spontaneous transition into a polarized, migrating state 
of a flat cell fragment. By changing some of the model parameters we can influence the 
degree of this symmetry-breaking instability, but so far we were not able to reproduce the 
observed longer-time stability of circular cell fragments, see again [63J. One reason for this 
failure seems to be, that in our model simplifications we assumed a constant distribution of 
freely diffusing myosin-II oligomers, contradicting experimental results on clearly expressed 
gradients of myosin-II concentration decaying towards the cell edge, see e.g. [60]. Therefore 
we have started to investigate the full coupled model system by including active tip poly- 
merization and the kinetics for myosin-II diffusion, binding and transport, in addition to 
the already implemented analogous kinetics and dynamics of integrin adhesion molecules. 
The next section presents the so far achieved results in the most simple but nevertheless 
quite instructive 1-dimensional situation. 

4.2 Induced onset of cell polarization and migration in the 1-D model 

We have simulated adhesion, polarization and migration of an idealized flat cell fragment 
having a fixed extension and only one degree of freedom to move, as can be observed 
experimentally, for example in Fig. [ijc). Thus, in a 1-D cross-section along its moving 
direction, the fragment is represented by an interval Q(t) = [xi,(t) , X{,(t) + L] of fixed length 
L, moving with body speed Vb(t) = Xb(t), so that the whole kinetics and dynamics within 
cytoplasm and dorsal membrane can be described by the corresponding 1-dimensional 
equations and conditions as in Section |3.2[ but now written in cell-centric coordinates, 
for example 0(t,y) = 6(t,Xb{t) + y) or v(t,y) = v(t,Xb(t) + y) — u&(t). Then, assuming 
the following particular restrictions, we have solved a simplified system of boundary value 
problems: 

1. The actin network is always sticky at the lamellar tips (no disruption), so that 6 > 
on the whole closed interval [0, L]. 

2. We assume active tip polymerization with a simultaneous parallel effect of the brown- 
ian ratchet {hq = 5 Pa/ fiM) and the clamp- motor mechanism {ttq = 5 Pa/ fiM). More- 
over, we suppose that the membrane-bound cortex shear viscosity strongly increases 
if myosin-II oligomers are bound: ^{mb) = 0.1(1 + 45m&) Pa min//im. For the actin- 
binding membrane proteins at the tip with maximal concentration A = 50 fiM and 
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Figure 7: Symmetric unpolarized state of 
a model cell (fragment) with 
length 10 [im. showing the con- 
centration profiles of 

(a) F-actin a, free and bound 
myosin-II m/ and nib, the force 
F onto the substrate and the ef- 
fective cortical F-actin velocity 
v c 

(b) of integrin proteins in the 
four different states, namely 
substrate- and actin-bound c sa , 
substrate-bound c s , freely diffus- 
ing Cf, and actin-bound c a . 



self-enhanced binding with dissociation constant K = 158.1 we use the pseudo- 
equilibrium 2a B (a) = A - ^(A - a) 2 + 2(A + a)K 2 /a + (K 2 /a) 2 + K 2 /a + a and 
a c = 0.1 a B . 

3. The dorsal membrane moves together with the cell (no slip at the tips), so that u = 
and v c = v = v + v\). The tip membrane tension difference is [tt]q = & u Lvb with the 
minimum always equal to a fixed positive constant tq = 25 Pa. Finally, the substrate 



force balance ( 39 ) reads 

-L 



<5> u Lv b = [ (f> c sa d(i + v b ). (42) 
J o 

This is an implicit equation to be solved for the migration speed v b , since the Neu- 



mann boundary conditions and the right hand side of the elliptic equation (19)) for 



v (after eliminating the pressure p) contain expressions that depend (linearly) on v b . 



Starting with the same unpolarized initial condition as in the previous Section |4.1[ we 
obtain a similar 1-dimensional symmetric configuration with no cell translocation and a 
central F-actin plateau, a central maximum of FA (focal adhesion) sites, i.e. c sa -integrins, 
as well as a centripetal F-actin flow, see Fig. [7j In addition, also the concentration of total 
myosin-II is enriched in the central region, with actin-bound myosin forming a central 
plateau, consistent with fluorescence pictures of unpolarized keratinocyte fragments, see 
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[63]. Moreover, the F-actin flow is highest at the boundary (with values V ~ 0.5 /im/min) 
due to the assumed active tip polymerization. 

In contrast to the previous 2-D model simulations, here we find robust parameter con- 
stellations yielding stability of this symmetric unpolarized state: even quite strong but still 
subthreshold perturbations of the F-actin concentration at one side induce only transient 
locomotion together with shifts in most concentration profiles: after some delayed over- 
shooting, the cell fragment returns to its non-moving non-polarized stable state in Fig. [7j 
However, if F-actin polymerization is continuously stimulated at one side (mimicking the 
effect of a chemo- or haptotactic gradient) then, not surprisingly, the cell fragment slowly 
polarizes and starts to persistently translocate in this direction (data not shown). 

Similarly, in order to mimic the mechanical stimulation experiments with keratinocyte 
fragments as performed in [63], we locally increase the activation (and F-actin-binding) of 
myosin-II at the left hand side for a certain time (up to 30 sec), which is thought to be 
analogous to the experimental push by a micropipette flow pulse, since by local compression 
of the cytoskeleton, filament alignment and thus myosin action will be enhanced. If time 
span or amplitude of this myosin-pulse stays subthreshold, the model cell responses only 
by a transient migration as described above and asymptotically returns to its stable resting 
state, see the speed curve in Fig. [8Fa). However, if the pulse strength exceeds a certain 
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Figure 8: Plots of cell migration velocity Vb over time after a perturbation of the unpolarized 
state in Fig. [7] by imposing an additional activation rate for actin-bound 
myosin-II oligomers over an induction period of 0.5 min locally at the left cell 
side: (a) with rate a+ = 20/min, showing a stable convergence of Vb towards the 
unpolarized speed zero, and (b) with rate = 50/min, showing the convergence 
towards a positive polarized speed ~ 0.012 /im/min. 

threshold, myosin-II and F-actin is condensed at the rear, whence active tip polymerization 
(V) is drastically reduced at this side, so that the cell rapidly starts to migrate in the other 
direction, see Fig.[8]jb). After cessation of the pulse, the migration speed is reduced, but the 
cell maintains its polarized locomotion state. This induced polarization, as a non-linear 
threshold behavior, is supported by a further positive feedback mechanism: Due to cell 
translocation, the focal adhesion sites (c sa ) are successively shifted rearwards relative to 
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the cell, which stabilizes the asymmetric polarization and later increases the locomotion 
speed to a constant asymptotic value (v;, ~ 0.12/zm/min). 




Figure 9: Migrating polarized state show- 
ing stable concentration profiles 

(a) of F-actin, a, free and bound 
myosin-II, mf and m&, the trans- 
duced locomotion force F and 
the effective cortical F-Actin ve- 
locity v c ; 

(b) of integrin-concentrations as 
in Fig. [7| Only FA integrins in 
the sa state are able to transduce 
force from the contractile actin 
network to the substrate. 

The migration speed is = 
0.125 ^m/min. 



In this migration state the cell (fragment) attains characteristic concentration profiles, 
see Fig. [9] Besides the already mentioned polar gradients of F-actin and myosin-II, the most 
impressive distribution is that of the FA sites: the c sa profile shows a characteristic broad 
peak behind the leading edge (as in the 2-D case above) , followed by a slight decrease and a 
second plateau of even more condensed adhesion sites in the back part of the cell (fragment) 
which, however, rapidly decays at the very rear, see Fig. (9^b). This last phenomenon is 
the theoretically expected and experimentally observed rear release of adhesion sites or 
integrins, not induced by any directed regulatory protein, but only by the fact that the 
rear part experiences a steep increase of the force F transduced to the substratum, see the 
plot in Fig. |9^a). By (37) this is proportional to the F-actin mass flow 0v c with respect to 
the substratum: While in the major front part of the migrating cell (fragment) the F-actin 
flow is retrograde and the centripetally pulling negative force is modest in amplitude, near 
the trailing edge the direction of flow reverses and the positive force becomes very strong, 
now centripetally pulling off the focal adhesion sites. Thus, the reason for cell translocation 
is not that there is "more adhesion" or "stronger force" at the front compared to the rear, 
as it is asserted in some models, cf. [17J. Indeed, for vanishing passive friction $ u due 



to (42) the total force integral even vanishes. The true physical reason for cell migration 



rather is the aforementioned asymmetry in the polarized cell state, expressed by a wide 
front region with modest rearward force and short rear region with strong forward force. 
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We finally investigate the dependence of migration speed on two cell physiologically 
important parameters, namely the adhesiveness of the substratum quantified by the relative 
number of available adhesion sites AdhO, and the responsiveness of the F-actin network 
measured, for instance, by the concentration ArpO of activated Arp2/3 proteins. These are 
mainly responsible for controlling F-actin polymerization by available free filament ends 
as they appear in the net assembly rate Q and the ratchet polymerization pressure (31). 
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Figure 10: Migration speed Vb of a model cell in its stable polarized state, plotted over 
varied parameters (a) of adhesiveness AdhO, relative amount of adhesion binding 
sites on the substratum (e.g. fibronectin coating) (b) of F-actin responsiveness 
expressed by ArpO [//M], the cytoplasmic concentration of activated Arp2/3 
complexes. 



The results depicted in Fig. 10 reveal the existence of optimal ranges for both parameters, 
consistent with earlier modeling results, see |17l I24j . and with experimental observations, 
see e.g. [271 [52]. 

Particularly, the migration response curve of Fig. [To^a) has been performed for fixed 
parameter ArpO = 10 fiM and for adhesion values < AdhO < 20, where according to 
Table [2] the adhesiveness proportionally influences not only the two adhesion rates a + = 
7 + = Adh0-a but also the passive membrane friction coefficient $ u = AdhO-6 Pa-min-//m -2 . 
This is based on the idealizing assumption that there are no other relevant exterior forces 
that would resist locomotion relative to the substratum, than those due to friction of the 
dorsal membrane with respect to an adhesive coating of e.g. fibronectin or collagen. 

Under these hypotheses and parameter choices our model simulations predict a minimal 
migration speed of Vb ~ 0.035 /xm/min for AdhO ~ 1, which increases towards the doubled 
speed when lowering the adhesiveness to almost zero. This surprising phenomenon is 
consistent with own experimental measurements of human keratinocyte polarization and 
migration |37j revealing slightly increased motility of polarized cells on glass compared 
to those on a low density fibronectin coat. The reason for this can be seen from the 



27 



corresponding profiles of FA concentration (c sa ) and substrate force distribution (F) for 
the two adhesiveness values AdhO = 0.01,0.2, see Fig. [IT|a,b). Due to ongoing retrograde 
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Figure 11: Substrate force distribution F and integrin concentrations in the four different 
states, plotted for varying adhesiveness parameter AdhO with corresponding 
migration speed as in Fig. [lO^a): (a) 0.01, speed 0.06^m/min, (b) 0.2, speed 
0.04^m/min, (c) 3.0, speed 0.125 /im/min, and (d) 20.0, speed 0.138 /im/min. 

flow of F-actin, working against internal viscosity and drag, the polarized cell (fragment) 
gathers the few FA sites on the back side in a way that the local FA and force distributions 
are almost symmetric (for almost vanishing adhesiveness AdhO = 0.01) but still with an 
additional negative (pulling) force plateau on the front side (though of tiny absolute value 
\F\ ~ 0.05 Pa). For increased but still small adhesiveness (AdhO = 0.2), the much more 
frequent FA sites start to accumulate at the rear, thus the asymmetry of forces is enhanced 
and the increased friction reduces the migration speed. 

On the other hand, for further increasing adhesiveness, AdhO > 1, the front plateau 
of pulling forces at the enriched FA 'carpet' is proportionally increased, see Fig. [Pi~|c,d). 
However, the dominant reason for the non-linear speeding-up response is the prominent 
increase of disruptive forces |F c (rear)| up to values of AdhO ~ 2, leading to a drastic 
reduction of resistive FA sites at the very rear, see the corresponding indicator curves 
in Fig. [l~2"|a). Only for very large adhesiveness, AdhO ~ 20, the FA sites again start to 
accumulate at the very rear in spite of a strong disruptive force there, as can be seen 
in Fig. [IT|d), causing a slight decrease in migration speed. Thus, while an adhesiveness 
between 5 and 12 leads to a saturated optimality in migration speed, within the range 
1 < AdhO < 4 the cell (fragment) has a high sensitivity for responding to an increase of 
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Figure 12: Variation of migration speed Vb and certain indicator functions at the cell edge 
in dependence of (a) adhesiveness AdhO and (b) F-actin responsiveness ArpO 
[/iM]: Plotted are the load \F C \ per FA site at the cell rear (bold) and tip (light), 
and the FA concentration c sa at the rear (bold, dashed) and tip (light, dashed). 
Ripples in the force curves reflect the slight stochastic noise that was added in 
the force equation (19). 
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adhesiveness with an up to 4-fold increase of speed. Carried over to the mean forward 
speed of competing leading lamellae in a whole cell model (as for the 2-D simulations in 
Fig. [6]), this could be used to explain polarization and haptotaxis of cells in spatial adhesion 
gradients, cf. [20J. 

In order to understand the analogous optimal migration performance for responsive- 
ness parameters in the range between 4//M and 10 fjM. of Arp2/3 concentration ArpO, see 
Fig. 10 'b), we plot the same indicator curves for FA concentrations and substrate forces in 
Fig. 12 b). Again, with increasing Arp2/3-induced actin polymerization there is an almost 



linear increase in pulling force at the tip (due to increased polymerization velocity V there) 
while the concave migration speed curve directly corresponds to the proportional curve for 
the disruptive forces at the rear and the resulting convex curve for the accumulated FA 
sites there. For larger ArpO values the growing accumulation of FA sites again leads to a 
slight speed reduction. Finally, the about 50% gain of speed for increasing responsiveness 
in the range 1 fjM. < ArpO < 4 /iM could again be carried over to a competing lamellar 
protrusion response in spatial chemotactic gradients, which are known to stimulate F-actin 
polymerization at the leading front. 

4.3 Migration speed in a simplified 1-D model 

In order to explore, which properties of our coupled cytoplasm- adhesion model are essential 
for obtaining the 1-D simulation results in the preceding section, we simplified the model 
by freezing the following variables and parameters, cf. 



1. Assuming fast diffusion and F-actin binding of free myosin-II oligomers, we take 
the pseudo-steady state condition ^ for bound myosin-II, but with the equivalent 
dissociation rate 5 m (a) = 5^ exp(— 2o/a op t) instead of ([8]), and obtain as contractility 
i/)(9) = tpo ■ 6 2 exp(— 2#/6> op t) with a coefficient ipo analogous to (41). 



2. Except the FA friction function $ we set all other friction coefficients (<p,$ u ) to zero. 



Thus, the cell migration speed is determined by the zero integral in (42). 



3. All 'exterior' pressures or tensions at the free membrane tip (7TQ,7TQ,Tr) are assumed 
to vanish, in particular, there occurs no active polymerization (V = 0). 

Starting with uniform FA concentration and a minor asymmetry in the distribution of 
actin, the system reaches a polarized and migrating state and the distributions of the 
various concentrations exhibit the same characteristics as in the full model (Fig [9]) . We 
again analyze the effect of substratum adhesiveness and F-actin responsiveness on the 



migration speed, see Fig. 13. In contrast to the observations for the full model (see Fig. 10), 
a saturation behavior of the migration velocity emerges with increasing adhesiveness instead 
of an optimal range. However, the responsiveness of the F-actin-network has an optimal 
value as in the full model. 
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Table 1: List of model variables and functions 
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Figure 13: Migration speed i>b for a simplified model cell in its stable polarized state in 
dependence of (a) adhesiveness or (b) F-actin responsiveness as described in 
Fig. 
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Further simulations after successive variation of the frozen parameters (see above) reveal, 
that an optimal range of adhesiveness could be gained only if a passive frictional load & u > 
was chosen to increase monotonically with adhesion strength AdhO, as it was assumed in 
the full model, see Tableland Fig. [To|a). Though such a friction appears to be physically 
realistic, the question remains whether such an optimal velocity response is of general 
biological relevance. For, the response of living systems is typically guided by adaptation: 
Migrating cells show the capability to dissolve their focal adhesions by proteolytic cleavage 
of integrins |55j which could be an adaptive strategy to optimize cellular adhesiveness 
on various substrates. On the other hand, the adhesiveness could be raised by active 
segregation of extracellular matrix proteins. 

Nevertheless, the universal optimality curve in dependence of F-actin responsiveness, 



seen in both Figs. 10 b) and 13 'b), seems to hold for any mechanical model, which imple- 
ments a saturating F-actin polymerization by association of regulating proteins as Arp2/3. 
This indicates that an optimal control of migration velocity by tuning certain chemical 
processes as F-actin polymerization, is more easy to be realized than tuning adhesion, 
and thus may have evolved as a generally effective strategy to regulate such mechanical 
processes as cell adhesion and locomotion. 



5 Discussion and outlook to further modeling 

Based on a larger set of interwoven mass and force balance equations for mechano-chemical 
processes within the cytoplasm and the surrounding plasma membrane, the presented con- 
tinuum model is a simplified though already sufficiently complex version of a more com- 
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Table 2: List of parameters for the 1-D simulations, if not otherwise specified. 
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prehensive 3-dimensional whole-cell model. From the 1-D and 2-D simulation results we 
conclude that even without organizing centers (as the cell nucleus or microtubuli, cf. |31j) 
or regulating systems (as the Rho/Rac control cascades, cf. [5l]) the F-actin cytoskeleton 
and its associated proteins having mechanical functions (myosin and integrin) constitute 
a self-organizing biophysical system with the ability of autonomous polarization and loco- 
motion. The transition form a symmetric, unpolarized stationary state into a polarized 
migrating state can be mechanically induced or spontaneous due to stochastic or chaotic 
fluctuations, depending on the (meta-)stability of the stationary state. Clearly, for repro- 
ducing the mentioned experimental data, the 2-dimensional system offers a wider spectrum 
of possibilities, one of which has been elaborated by Kozlov and Mogilner [33J. They prove 
a defined bi-stability between the radially symmetric cell shape and a polarized, circu- 
larly indented shape by allowing for different anisotropic organization of the actin-myosin 
cytoskeleton. In comparison to this, further numerical 'experiments' using our coupled vis- 
cous flow-transport-diffusion-reaction system with free boundary should be performed to 
explore the capability of isotropic continuum models. Though local orientation or alignment 
of actin filaments will have an enforcing effect onto the dynamics, interactions and feedback 
between F-actin, cross-linkers, motor-proteins and adhesion complexes (for a 1-dimensional 
model see [21])) already the presented simulation results show that internal concentration 
gradients and directional flow, leading to persistent polarization, can emerge from local fluc- 
tuations. Thereby super-threshold pattern formation arises, not as in common excitable 
media by a purely reaction-diffusion feedback, rather by a spatially distributed coupling 
between chemically induced force generation/relaxation and mechanically induced bond 
formation /dissociation. 

More than one decade ago, Bereiter-Hahn and Liiers [8] claimed that, based on their ex- 
periments and measurements with migrating keratocytes, polarization and directionality of 
cells are determined by local variations in actin network stiffness and hydrostatic pressure. 
Together with the more recent high-resolution measurements of dynamic vector field (as 
F-actin velocity, forces onto pliable substrates, or directionality of FA shape deformations), 
it could soon be possible to probe and evaluate the various assumptions and hypotheses 



on mass and force balance conditions, which we postulate in Sections 3.1 and |3.2| For 
example, the biphasic correlations between retrograde actin flow speed and traction stress 
detected by Gardel et al. [22] could well be resolved by the measured and simulated FA 
gradient away from an advancing lamellar tip, and by our central modeling assumption 
that the local traction force is proportional to the product c sa ■ v c of FA concentration 
and cortical F-actin velocity. 

Moreover, the distribution of flow and force, extensively discussed for the 1-dimensional 
case (in Section |4.2|), has to be quantitatively characterized for the 2-dimensional free 



boundary situation (Section 4.1), where the overall picture is the same: Wider distribu- 
tion of traction forces in the leading front region (with maximum flow speed at the tip, 
depending on the active polymerization pressure) and strong disruptive counter-forces in 
a confined region of 'rear release'. The actin flow patterns of our 2-D simulations show in 
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principle the same characteristics as experimental data by Yam et al. |64J for fish kerato- 
cytes. During the process of polarization, the simulated cell looses its radial symmetry of a 
uniform inward flow pattern, leading to a stronger inward flow in the rear region, whereby 
the most pronounced forces come from the two sides relative to the establishing locomotion 
direction, see Fig. [6^c). Indeed, traction force experiments on polarized fish keratocytes 
have revealed the same pattern, with major traction forces from the two flanks at the rear 
and weaker retrograde forces at the front [50!. 

Further evaluation of the presented model should compare the simulated integrin concen- 
tration profiles with more precise quantitative data on the temporal turnover and spatial 
distribution of FAs in migrating cells, as already mentioned above. For this reason, the 
model has to be extended to include stochastic clustering dynamics of integrins, which 
probably is an essential feature of FAs determining the detailed patterns of F-actin flow 
and traction forces. We emphasize, that our model results do not depend on the specific 
mechanisms of myosin-II contraction. Essential for the emergence of coupled flow and mass 
gradients is the bimodal (cubic- like) stress function S = S(6) in Fig. [5] which similarly ap- 
pears in a modified mechanical situations as, for example, the nematode sperm motility 
system, see [D]. 

The experimental response curves for the migration speed of polarized cells depending 
on adhesiveness [52] show a broad optimum range with a logarithmic slow decay for larger 
AdhO values and a steep decay for lower ones, similar to the curves found in our 1-D model, 
Fig. 10 and in earlier analogous simulation results using alternative visco-elastic mechanics 
|17} [23] . The deviating behavior for very low adhesiveness values hat still to be investigated 
with a more precise quantification of other (weaker) substrate-mediated forces that may 
come into account. 

An optimum curve for the protrusion speed in dependence of the barbed end density 
B, thus also of the concentration ArpO, has been obtained by Mogilner and Edelstein- 
Keshet |42j with the aid of 1-dimensional diffusion-reaction-transport systems, including 
membrane resistance and the brownian ratchet mechanism. However, they do not model 
retrograde actin flow or force transduction to the substratum, so that the shape of their 
curve ([12], Fig. 6) essentially differs from our results in Figs. 10 Jo) and 12 b). 



Finally, in generalization of the model simulations by Stephanou et al. [59] and in compar- 
ison to other modeling approaches mentioned before, our full 2-dimensional model should 
be used for reproducing the experimentally observed translocation paths and deformation 
structures of single blood or tissue cells in culture, which show well-expressed phases of 
persistent polarization and locomotion, interrupted by events of speed reduction, contrac- 
tile rounding, successive re-polarization and migration. Statistical analysis, in the same 
spirit as having been performed for the original annual lamella model [3], could help to 
reveal the microscopic mechanisms that are responsible for cell motility behavior on the 
macroscopic level. 
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